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Polar molecules, in strong-field seeking states, can be transported and focused by an alternating 
sequence of electric field gradients that focus in one transverse direction while defocusing in the other. 
We show, by calculation and numerical simulation, how one may greatly improve the alternating 
gradient transport and focusing of molecules. We use a new optimized multipole lens design, a 
FODO-lattice beam transport line, and lenses to match the beam transport line to the beam source 
and to the final focus. 

We derive analytic expressions for the potentials, fields, and gradients that may be used to design 
these lenses. We describe a simple lens optimization procedure and derive the equations of motion 
for tracking molecules through a beam transport line. As an example, we model a straight beamline 
that transports a 560 m/s jet-source beam of methyl fluoridel5 m from its source and focuses it to 2 
mm diameter. We calculate the beam transport line acceptance and beam survival, for a beam with 
a velocity spread, and estimate the transmitted intensity for specified source conditions. Possible 
applications are discussed. 

PACS numbers: 39.10.+j, 33.15.Kr, 07.77.Gx 



X 



I. INTRODUCTION 

A polar molecule has an intrinsic separation of charge 
on which an electric field gradient exerts a force. The 
force, F x , in the (transverse) x direction is: 



F, 



dW 

dx 



dW dE 
~dE~dx 



(1) 



where W is the potential energy of the molecule in an 
electric field E (Stark effect) of magnitude E = (E 2 + 
Ey) 1 / 2 ; and similarly for the force, F y in the y direc- 
tion. The force and gradient are in opposite directions 
(weak-field seeking) for rotational states that become less 
tightly bound in an electric field (dW/dE > 0), while the 
force and gradient are in the same direction (strong-field 
seeking) for rotational states that become more tightly 
bound in an electric field (dW/ dE < 0) . The J = state 
is always strong-field seeking and all rotational states be- 
come strong-field seeking in the limit of strong electric 
field as shown in Fig. [l]. 

Focusing a beam of molecules, traveling in the z di- 
rection, is achieved using static two-dimensional (x, y) 
electric field gradients. We neglect end field effects and 
assume E z = inside the focusing lenses. Molecules 
in weak-field seeking states can be focused, in both di- 
rections, using quadrupole and/or sextupole fields that 
have a minimum in the electric field in both directions. 
Molecules in strong-field seeking states, however, can be 
focused in only one transverse direction, while being de- 
focused in the other, because it is not possible to have a 
maximum in the electric field in both dimensions (in free 
space). 



Molecules in strong-field seeking states have been 
transported and focused by an alternating sequence of 
electric field gradient lenses ]l|, ||, ||, f§ ||, ||, @] (as have 
neutral atoms j8|), but it has been neither as success- 
ful, nor as widely used, as has quadrupole and sextupole 
focusing for molecules in weak-field seeking states ^ ^| . 

In this paper we show, by calculation and numerical 
simulation, how one may greatly improve the alternat- 
ing gradient transport and focusing of molecules by op- 
timizing the lens field geometries so that nonlinearity in 
dW/dE is compensated by dE/dx to produce a linear 
restoring force. We use this optimized multipole lens de- 
sign in a FODO-lattice beam transport line and use ad- 
ditional lenses to match the beam transport line to the 
beam source and to the final focus. 

The remainder of this paper is organized as follows: 
Section || discusses alternating gradient focusing, trans- 
port lattices, matching lenses, and linear optics. Section 
[II derives the formulae for designing linear focusing ele- 
ments and presents examples of lenses. Section |y] derives 
the equations for molecular beam transport and models 
both a simple 30 m-long FODO lattice and a complete 
15 m-long trans port l ine with entrance and exit matching 
lenses. Section [IV G| estimates the intensity of a methyl 
fluoride jet-source beam, transported through a beamline 
and focused. Section ^ discusses the use of strong-field 
seeking states, and possible applications of alternating 
grading focusing and transport through very long beam- 
lines. 



II. ALTERNATING GRADIENT BEAM 
TRANSPORT AND FOCUSING 
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A beam of charged particles can be focused and trans- 
ported over almost unlimited distances by alternating F 
and D type magnetic quadrupole lenses. The F type lens 
focuses the beam in the horizontal (x) direction while 
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FIG. 1: Energy levels of low-lying rotational states of a 
diatomic or symmetric-top molecule in an electric field (Stark 
effect). The levels were calculated using the rigid rotor model. 
The horizontal axis is in dimensionless units of ui = d e E/B, 
where d e is the electric dipole moment, E the electric field, 
and B the rotational constant. The vertical axis is in units 
of B. Energies for small ui are shown in (a) and energies for 
large u in (b). The potential energy of the molecule, W, is 
the change in energy from zero electric field and for very large 
ujW — » d e E. The levels that descend are strong-field seeking 
(see text). For the J = state of methyl fluoride (CH3F), in 
the K = state, w = 1 at 2.78 MV/m. For the J = state 
of CsF, lo = f at 0.135 MV/m. 



defocusing it in the vertical (y) direction. The D type 
defocuses the beam in the horizontal direction while fo- 
cusing it in the vertical direction. 

A complete alternating gradient transport line begins 
with a beam source whose output is optically matched, by 
lenses, into the acceptance of a transport section which 
in turn is matched, by lenses, to a final focus. The final 
focus can be, for example, at an experimental target or 
at a matching point for injection into a (storage) ring 
lattice. A complete transport line for molecules, starting 
with a jet source and skimmer, is shown in Fig. 0. For 
long-distance beam transport, the lens system providing 
the highest beam acceptance and requiring the lowest 
focusing strength is a lattice of identical FODO-cells, in 



which the F/D- lenses are separated by drift spaces (O). 

To match the molecular beam source to the transport 
line, and to match the transport line to the final focus, we 
use doublet (FD/DF) and triplet (FDF/DFD) type lens 
configurations. These lens configurations generally re- 
quire higher focusing strengths than do the FODO cells. 
Doublets are typically used for asymmetric focal points 
(unequal x and y dimensions) while triplets are used for 
symmetric ones. In the F/D doublet, a net focusing in x 
and y occurs because particles first focused in x are, in 
the second element, closer to the axis and therefore less 
defocused. Particles first defocused in y are, in the second 
element, further from the axis and thus more strongly fo- 
cused. 

In all lenses, linear focusing is needed for optimum 
optical properties. Linear focusing requires that the 
force in each direction be linearly proportional to the 
displacement in that direction (|_F x /a;| = constant and 
\F y /y\ — constant), and independent of the displacement 
in the other direction (uncoupled motion). Strong non- 
linearities in focusing elements will result in loss of beam, 
generation of beam halo, and growth in the transverse 
emittance (the product of the angular divergence and 
the spatial dimension) producing larger beam sizes. 



III. LINEAR OPTICS 

A. Potential energy of a molecule in an electric 
field 

To determine the lens shape that will produce the 
most linear force on the molecule, we need dW/dE for 
Eq. (|l|). This quantity will change with the J, \mj\ rota- 
tional state and the electric field strength. If a number 
of rotational levels have a similar dW/dE, then one lens 
design will be nearly optimum for all of them. In the 
limiting case of a strong electric field (large to in Fig. [I] ), 
dW/dE = -d e . 

For polar molecules in weaker fields, we calculate the 
interaction energy, in the rigid rotor approximation, fol- 
lowing the approach of von Meyenn JTo| . The Hamilto- 
nian operator is H — BJ 2 — d e E cos 6, where B is the 
rotational constant, and the direction cosine matrix ele- 
ments, which couple (J, mj) with (J+l, mj) and (J— 1, 
toj), are taken from Townes and Schawlow We di- 
agonalize the matrix for terms through J = 30. The first 
few levels are shown in Fig. |l|. As u> = d e E/B and B 
determine the Stark effect for each J, \mj\, it is straight 
forward to construct a simple function for the Stark ef- 
fect for any J, \mj\. For the J = state of a diatomic 
molecule or symmetric top (K = 0) molecule, a satisfac- 
tory approximation [JL2| is: 



W(E) = 



duj 2 B 



C x d\E 2 



1 + C 2 uj B + C 2 d e E 



(2) 



with d = -0.2085 and C 2 = 0.2445. This expression 
works best for small and intermediate values of u>. 
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FIG. 2: Schematic diagram of an alternating gradient beam transport line for molecules, 
skimmer, beam matching optics, FODO lattice, and final focus optics. 



It consists of a jet source and 



From Eq. (||), for the J = state, the field derivative 
of the potential is: 



dW 
~dE 



W 
~E 



2B + C 2 d e E 



B + C 2 d e E 



(3) 



Similar expressions may be found for other rotational lev- 
els. A perturbation expression for W can be used for ro- 
tational states in (very) weak electric fields [[ll] . In some 
applications, using dW/dE from either the strong field 



limit or the weak field limit will be sufficient to design a 
linear lens. 

Finally, we note that all atoms and molecules, includ- 
ing nonpolar molecules, are polarizable with an interac- 
tion energy W a — — \aE 2 , where a is the dipole po- 
larizability For laboratory electric fields this 

interaction is much smaller than the interaction with a 
molecular electric dipole moment, but it can be used to 
focus atoms [pi or decelerate them [15l. 



J 



B. Electric field gradient of a focusing lens 



Any lens electrostatic potential, <& can be expressed as the following multipole expansion in cylindrical coordinates, 
r, 9: 



9) = E r 



°° h / r \ n °° n / r\ n 

co S (n9) + J2-(-) sin(r 
, n \r Q J ^ n \r J 



where Eq is the central field [E = E$(b\ + af) 1 / 2 ], tq is a scaling length, and b n (a n ) are the dimensionless constants 
of the 2n- pole strengths for normal (skew) elements. 

Setting b n = a n for simplicity, converting to Cartesian coordinates, and retaining only the oi, &3, and 05 terms (see 
section 



III C[ ), the lens potentials for normal multipoles are: 

a 3 (x 2 ~3y 2 ) a 5 (x 4 - 10x 2 y 2 + 5y 4 ) 



$ N (x,y) = E x 



a 1 



3rg 



5r 4 



or for skew multipoles are : 



a 3 {3x 2 ~y 2 ) a 5 (5x 4 - 10x 2 y 2 + y 4 ) 1 



3r 2 



5r 4 



where Eq is the field on axis for aj. = 1. 

Both of these potentials give the same total electric field (E = — V$) 



(4) 



(5) 



E(x,y) = (E 2 x +E 2 )^ 2 =E 



2 2a ia3 (x 2 - y ) a 2 (x 4 + 2x 2 y 2 + y 4 ) 2 ai a 5 (x 4 - 6x 2 y 2 + y 4 ) 

a l H 15 1 14 1 14 



1/2 



(6) 



which has the electric field gradients: 
dE _ 2a 3 E 2 G x x 



dx 
dE _ 
dy 



rl E 
2a 3 E 2 G y y 
rl E 



where 



G x (x,y) =ai + 



(7) G y (x,y) = ai- 



a 3 {x 2 +y 2 ) 2 ai a 5 (x 2 - 3y 2 ) 



a 3 r 2 



a 3 {x 2 +y 2 ) 2 ai a 5 (3x 2 - y 2 ) 



a 3 rl 



(8) 
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FIG. 3: Equipotentials (truncated) in x and y of skew F 
lenses. The actual electrode can be fabricated from any set 
of equipotentials and, if dimensions are in mm, they match 
lenses used in the examples in Section [f^. A two-wire field 
lens (with multipole components: a± = 1, az = —1, as = 
1 . . .) is shown in (a). An optimized multipole lens (with 
multipole strengths a\ = 1, 0,3 = — 1, as = 0) is shown in (b). 
The lens in (b) is designed to focus in the strong field limit. 



C. Force on the molecule due to an electric field 
gradient 

The force on a polar molecule in the x or y direction 
can now be calculated using Eq's. (0, |?], and |s|) with a 



suitable expression for dW/dE from section [II A . 

Non-linear forces, inside the focusing lenses, limit the 
maximum beam size that one can transport without suf- 
fering beam losses, emittance (size) growth, and beam 
halo. Non-linear forces generally arise from higher-order 
multipole components of the lens electrostatic potential. 
Cylindrical electrodes, which are two-wire field lenses 
(with ro = half-gap), shown in Fig. || (a), have long 
been used to focus molecules in strong-field seeking states 
||, |, |, |, §|, . They contain multipoles of all odd or- 
ders and of equal strengths ( a% = 1, as = —1, 05 = 1 . . .). 
As we will see below, its strong decapole field (ag = 1) 
reduces the area of the lens over which the focusing is 
linear (dynamic aperture). A beam transport line using 
two-wire field lenses will be limited to a smaller diameter 
beam compared to a similar transport line using opti- 
mized multipole lens such as the one shown in Fig. ^ (b). 

To produce a larger linear region we try to limit the 
multipole fields to those that are essential for the optics. 
We include a dipole field (ai) to give a non-zero field on 
axis, without which strong-field seeking molecules defo- 
cus in both transverse directions; a sextupole field (03) 
to provide the linear focusing force (of order r 1 ); and 
a decapole field (05) to correct the non-linear forces (of 
order r 3 ) produced by the sextupole field. We omit the 
quadrupole field (02), because it bends the beam, and 
we omit the octupole field (04), because it introduces 
stronger non- linear forces (of order r 2 ) than those of the 
sextupole/decapole fields. All other multipole strengths 
are set to zero. In a real lens, however, the electrodes 
may be truncated equipotentials, with the consequence 




FIG. 4: Contour plots of the force in the x direction (F x ), 
on a molecule inside a lens, as a function of x and y position. 
The lines are contours of equal F x in arbitrary units. In this 
example we use the large Stark effect limit (W = —d e E) and 
ro = 6 mm. (a) shows F x for a two-wire field lens, (b) shows 
F x , for an optimized multipole lens (as = 0). A lens that was 
completely linear in F x would have vertical contour lines with 
uniform spacing. 



that small residual higher order multipoles will remain. 

Since the lens potentials, as defined by Eqs. (||, ||), have 
two free parameters, we chose a\ = 1, making Eq the 
central field and we choose 1 0.3 1 = 1, making it easier to 
compare the optimized multipole lens to a two-wire field 
lens. The remaining choice is the decapole field strength 
as, which we use to optimize the lens optics. We do this 
in two ways: first by calculating, for specific molecules, 
the forces F x and F y inside a lens and comparing their 
linearity in x and y for different 05 , and second, in section 
IV, by simulation of beams in model transport lines. 



If we calculate the horizontal force, F x on a molecule 
whose potential energy \sW — —d e E (strong field limit), 
we find that it is most linear on-axis (y = 0) for no de- 
capole field (05 = 0). For this example, the constant- 
force contours for F x in the x, y plane are shown, for the 
two wire field lens, in Fig. [l] (a) and, for the zero-decapole 
lens, in Fig. |(b). The zero decapole lens is seen to have 
far better linearity than the two-wire field lens. (This is 
also true for F y which is not shown here.) In a similar 
way, we find that the optimized multipole lens outper- 
forms the two-wire field lens in the low electric field limit 
(quadratic Stark effect), where the most linear field on 
axis has 05 = — A. 
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D. End effects 



reduces to: 



Real lenses have finite length, and the two-dimensional 
fields inside become three-dimensional fields at the ends. 
There are three effects from this. First, the z-component 
of the field gradient affects the transverse motion. This 
will be small if the lens aperture is small compared to 
the physical length of the lens. 

Second, the lens field acts over a distance that is dif- 
ferent than the physical length of the electrodes because 
of (extending) fringe fields at the ends. For linear design 
optics we approximate the real lens by a lens of constant 
central field Eq over an effective length, different than 
the physical length of the lens. The effective length is 
determined by integrating the actual central field (found 
by numerical calculation or measurement) through the 
lens, including the ends, and dividing by Eq. 

Third, the z-component of the electric field gradient, 
present at the ends of the lens, has an effect on the z- 
component ( longitudinal) velocity. The non-relativistic 
Hamiltonian (neglecting gravity) for a molecule of mass 
to and total velocity v is: H = \vnv 2 + W. The Hamil- 
tonian is conserved if the electric field is static. Conse- 
quently, a molecule in a strong-field seeking (weak-field 
seeking) state will gain (lose) kinetic energy entering the 
lens and then lose (gain) the same amount of kinetic en- 
ergy upon exit [jl5|, lf|. For the central (x = y = 0) tra- 
jectory, the longitudinal velocity, v z e, of the molecule, in 
the lens is: 



x + Kqx = y" — Kay = 
where Kq is the lens linear focusing strength: 



(10) 



Kq 



2a 3 



Bo 
— j- 



dW 



dE 



2a 3 rfe-Eo 
~>T mv 'zE 

a 3 aEl 

1 1 



djK Ci(2B+C 2 deE ) 

mV l E 



for W 



-d P E 



1 9 

-2 aE 



(B+C 2 d c E ) 



for W a = 
for J = 0, from Eq(EL) 



For a molecule in a strong-field seeking state 
(dW/dE < 0), the lens will focus in the x direction (F- 
lens) for 03 < 0, and focus in the y direction (D-lens) 
for (13 > 0. For a molecule in a weak-field seeking state 
(dW/dE > 0) the F and D directions are reversed. The 
lenses can also be used to focus atoms (dW a /dE = —aE) 
which, in their ground states, are always strong-field 
seeking. 

The linear design optics will determine the lens focus- 
ing strength Kq. Then, the choice of scaling radius ro, 
will give the required lens field, Eq, from Eq. (|IT|). Small 
ro allows us to use low electric fields but, as seen from 
Eq's. (||, ||), large ro is needed for increased linearity. The 
value of tq chosen will then depend on the maximum elec- 
tric field strength and the beam size to be transported. 



v 2 zE 



W/2m 



where v z q is the velocity of the molecule in the drift 
spaces. For most thermal or jet-source molecules in lab- 
oratory electric fields, the change in velocity is a small 
effect. (A beam of 560 m/s methyl fluoride molecules in 
the J = state increases its velocity by 0.2 m/s upon 
entering an electric field of 10 MV/m.) 



IV. BEAM TRANSPORT 

A. Equations of motion 

To track the trajectories of molecules passing through 
combinations of lenses, we use the transverse nonlinear 
equations of motion for a molecule traveling in the z di- 
rection, given from Eq's. (|l|, |7|, and ||): 



2a, 



El G x x dW 



'0 " LU ZE 
?2 



y 



2a 3 El G y y dW 



J ZE 



E dE 



(9) 



with trajectory gradients defined as x' — dx/dz — 
v x /v ZE and y' = dy/dz = v y /v Z E- 

We start the design optics using completely linear 
lenses (G x = G y = 1, and E = E ), for which Eq. (^) 



B. Example of methyl fluoride 

For a realistic beam transport simulation, we assume a 
beam of methyl fluoride (CH3F), in the J ~ rotational 
state, having a longitudinal velocity vzq — 560 m/s. 
(This is the approximate velocity of a beam produced 
by seeding methyl fluoride in an argon jet s ource with a 
reservoir temperature of 300 K (See section IV G[ ) . The 
electric field derivative of the potential energy (of the J — 
state), is given by Eq. (|) with d e = 6.25 x 10" 30 J/V/m 
(1.86 Debye) and rotational constant B = 1.76 x 10~ 23 J 
(0.88 cm' 1 ). 

To study long distance transport, we model a FODO 
lattice consisting of a sequence of identical F and D lenses 
separated by drift spaces (O). Then, to study a complete 
beam line, we add an upstream lens section, to match 
the beam from the jet source into the FODO lattice, and 
a downstream section for a final focus of the beam (see 
Fig- 1). 



C. FODO lattice 

We chose a simple FODO lattice consisting of identical 
FODO cells (see Section ||). The optics of a FODO cell 
starts in the center of an F lens (or D lens) and ends at 
the center of the next F lens (or D lens). At the ends of 
the FODO cell the beam is at a waist (defined as zero 
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slope in the beam size) with a maximum size in one plane 
and a minimum size in the other plane. 

For this study, we use 10 cm-long F/D lenses, with scal- 
ing length tq — 6 mm, separated by 40 cm drift spaces, 
giving a FODO cell length of 100 cm. This design leaves 
80 % of the cell unfilled. It is economical to build, but 
does not have as large an acceptance as designs that filled 
more of the lattice (see section [v]) . 

In a FODO lattice the motion of a particle is peri- 
odic in phase-space. The phase-advance measures how 
far along the period it has proceeded from its initial start- 
ing point. The transverse linear optics are characterized 
by the phase advance in the FODO cell /i c . For our cell 
we choose /i c = 7r/3 rad, for which all particles return 
to their initial phase-space position after 6 FODO cells. 
This specifies a central field of Eq — 3.23 MV/m, which 
is not close to breakdown, as well beam size that is 
nearly the minimum possible. 



D. Modeling beam transmission 

1. beam distribution in phase space 

We take the molecular beam to be continuous (un- 
bundled) and monoenergetic, in which case, it is com- 
pletely specified by its density in the [x, x', y, y') phase 
space (distribution function). The beam size in a trans- 
port line depends, not on the density, but on the four- 
dimensional volume (V4,d) occupied by the beam. This 
can be defined, experimentally, in a number of ways: (a) 
the volume containing some fraction of the beam, (b) 
the volume calculated from the root-mean-square (rms) 
beam sizes, or (c) the volume defined by a set of col- 
limators. These four dimensional volumes can be re- 
lated to an equivalent theoretical distribution function 
of Kapchinsky and Vladiminsky (KV) [[l7). With the 
KV distribution function, we can calculate the linear fo- 
cusing beam sizes (a x ,a y ) along the beam transport line 
and, in so doing, model the design optics. 

The KV distribution consists of a uniform density of 
particles on a hyper-ellipsoid in four-dimensional phase 
space. At a beam waist (a x = a' y = 0), the hyper- 
ellipsoid is given by: 



K a a y = 



t KVy 



(12) 



and has the volume, Vkv = \it 2 a x a x i a y a y i . Typically 
we have beam waists at the beam source (minimum), 
at any focus (minimum), and inside the focusing lenses 
(minimum or maximum). 

The beam sizes along the beam transport line are given 
by the uncoupled envelope equations, which for the equa- 
tions of motion [Eq. (|lC|)] are: 



K a x 



,2 

e KVx 



where exvx and tKVy are the invariant transverse emit- 
tances which, at a waist [see Eq. (|l^)], are simply given 
by e-KVx = a x a x f and e KVy = a y a y > . 

If we project the KV distribution onto the (x, x') plane, 
we obtain a uniform density of particles inside an ellipse 
of constant area irtKvx- For an arbitrary beam distri- 
bution, the equations of motion, for the rms beam sizes, 
have the same form as Eq. (|lj). This defines the rms- 
equivalent KV distribution sizes as a x = 2cr x , a x i = 2a x >, 
a y = 2(7 y , and a y > — 2oy . Under linear forces, a molecule 
remains on the same KV surface on which it started, with 
the shape of the ellipse changing but its volume remain- 
ing constant. 

Since the non-linear forces inside a lens become 
stronger, the further the molecule is from the central axis, 
lenses may be evaluated by computing the increase in 
non-linear effects with increasing beam size. For this we 
use nested KV distributions of increasing volume, each of 
which is characterized by its volume, Vkv- The simula- 
tion results are then independent of the characteristics of 
the initial beam and can be applied to non-surface beam 
distributions. 



2. beam survival 

We calculate the survival of a 560 m/s methyl fluoride 
beam in the J = state, through a 30-m FODO lattice, 
as a function of the decap ole strength (05). The lattice is 

and we assume that the beam 



described in section [V C 



is already matched to the lattice. As molecular beam jet 
sources are usually axi-symmetric, we take the emittance, 
ckv to be the same in both transverse planes. Then, the 
initial matched beam sizes and divergences, for the 7r/3 
phase-advance FODO cell, are: 



Q>max &X 

a x i 
a y 

Oy> 



{^KvPmax) 1 ^ 2 
{^KvPmax) 1 ^ 2 



a/2 



where m i n = 0.587 m and f3 max = 1.703 m. As is in- 
dependent of the emittance, we use the initial maximum 
beam size, a max as our KV distribution size parameter : 



1 



i - -2,2 " "-max 

Vkv = 2^ € kv = 



2/3. 



2 

■max 



The relation between the beam size and the emittance is 
shown in Fig. %. 

The trajectories were numerically integrated along the 
FODO lattice using a 4'^-order Runge-Kutta algorithm. 
We take as lost those particles whose transverse displace- 
ment, in the beam transport line, in x and/or y became 
too large. The transported particles, which started on a 
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maximum envelope size (mm) 



FIG. 5: Emittance versus beam envelope size of the shell in 
the Kapchinsky-Vladiminsky (KV) equilibrium beam distri- 
bution. 



zero-thickness KV-surface, finished up on a smeared-out 
fuzzy one, producing emittance growth and halo. 

From the simulations, we found the fraction of the ini- 
tial methyl fluoride beam surviving as a function of po- 
sition along the beamline. We did this for two different 
lens designs (two- wire field lens and optimized multipolc) 
and for selected initial KV beam sizes. The results are 
shown in Figs. || (a) and ^| (b). 

For the two- wire field lens, the largest KV beam that 
can be transmitted without loss is 1.4 mm. For the op- 
timized multipole lens, the largest lossless KV beam is 
2.5 mm. If the initial phase-space density of the beams 
is approximately constant, the relative beam intensity is 
given by the ratio of the lossless phase- space volumes. 
This is a factor of (2.5/1.5) 4 rs 10 improvement by using 
the optimized multipole lenses in place of the two-wire 
field lenses in this FODO lattice. 

Figures H (a) and|| (b) also show that most of the beam 
losses occur in the first 5 m and by 30 m the losses are 
essentially complete. This suggests that extending the 
beam transport line to much longer distances will not 
further reduce beam survival. 

To find the most appropriate values of the lens de- 
capole strength, for other molecules, we also studied 
beam survival for the two limiting cases of strong elec- 
tric field (linear Stark effect) and weak electric field 
(quadratic Stark effect). A KV beam size of a max = 2.5 
mm was used and the results, along with the results for 
methyl fluoride, are shown in Fig. |7|. 

For the strong field limit, as well as for the methyl 
fluoride example, we obtain the highest transmission for 
as = 0. In the weak electric field limit the optimum value 
of 05 is shifted to 05 ~ —0.2, this nonlinearity in dE/dx 
compensating the nonlinearity in dW/dE. Thus, an a$ 
near zero would be a good choice for CsF molecules in 
low-lying rotational states, and an 05 near -0.2 would be 
a good choice for ground-state atoms. 
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FIG. 6: Calculated beam survival for 560 m/s methyl fluoride 
in the J — state in a 30 m-long FODO lattice constructed 
of (a) two wire lenses or (b) optimized multipole lens with 
as = 0. In each case ro = 6 mm and the central field, Eq — 
3.2 MV/m. Survival is plotted as a function of (KV) beam 
size. For the 100% beam survival, the beam intensity scales as 
KV beam size to the fourth power (see text). Consequently, 
the FODO lattice of optimized multipole lenses will transport 
about a factor of ten more beam than the same FODO lattice 
of two- wire field lenses. 



E. Velocity dependence of the beam transmission 

So far we have considered only a monoenergetic beam. 
To look at the tolerance of the different lenses to en- 
ergy deviations, we calculated the transmission for ini- 
tially identical KV distributions but with different ener- 
gies (560 m/s being the matched velocity.) 

The results are shown in Fig. ||, for KV beam size of 
o-max — 1.5 mm, for both the two- wire field lens and 
for the optimized multipole lens. Both do quite well. At 
energies from 0.65 to 1.2 times the nominal energy of 640 
K (velocity 560 m/s), transmission in the FODO lattice, 
of two-wire field lenses, is 50% or more of its maximum 
value. The FODO lattice of optimized multipole lenses 
does even better, as we would expect, since it has lossless 
transmission up to a KV size of a max — 2.5. One should 
note that Fig. g does not represent the energy acceptance 
of a complete transport line, which will be limited by the 
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FIG. 7: Survival of 560 m/s monoenergetic beams through 
the 30 m FODO lattice for a KV beam with a ma x = 2.5mm 
as a function of decapole constant (as). Results for a methyl 
fluoride at 3.2 MV/m are shown as open squares, the limiting 
cases of very strong and very weak electric field as open circles 
and filled diamonds respectively. 
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FIG. 8: Calculated transmission, as a function of relative 
beam kinetic energy, of a beam of methyl fluoride in the J = 
state through a 30-m FODO lattice. The nominal energy is 
640 K (560 m/s) and the KV beam size is a ma x = 1.5 mm. 
Points for the FODO lattice using two-wire field lenses and 
optimized multipole lenses (with decapole constant as = 0) 
are shown as squares and circles, respectively. 



source and matching optics. 



F. Beam matching and point focusing 

We complete the model transport line by adding 
matching optics upstream and downstream of the FODO 
lattice. This matches the source to the lattice and the 
lattice to the final beam focus. The matching, in both 
cases, is achieved using a doublet- and a triplet- lens con- 
figuration as shown in Fig. |[ As the matching lenses re- 
quire stronger focusing, producing more non-linearities, 
we raise the lens scaling length in the matching lenses to 
tq = 12 mm. And since most of the beam losses occur 



early in the FODO lattice, we shorten the lattice to 15 m 
to reduce the computing time necessary for the simula- 
tion. The beam envelope, in x and y, and the placement 
of the lenses, is shown in Fig. |^, for an initial beam which 
is close to the linear focusing limit in the lenses. 

For our example, of a 560 m/s beam of methyl fluoride 
in the J = state, we assume that the jet source has a 
very small orifice and a skimmer, of 1.5 mm diameter, 
is placed 25 cm upstream of the first lens. An angular 
spread of ± 3.3 mrad is fixed by collimators (see Fig. ||). 
This gives an initial beam emittance of 2.5 mm-mrad in 
both transverse planes, and a maximum beam size of 2.1 
mm in the FODO lattice. In addition, the initial beam is 
taken to have a Gaussian energy spread of ± 10 % (rms). 
The calculated beam transport, through this beamline, 
is 84% of the entering beam, the losses being due to the 
energy spread. Thus, we can transmit and focus most of 
the velocity distribution from the jet source. (If needed, 
the beam can be focused, at the end, to an even smaller 
size. However in this case, we have to increase the size 
of the beam in the doublet/triplet lenses, where the non- 
linear forces will produce emittance growth and beam 
halo. If we arc not limited by the electric field strength, 
we can correct this by using lenses with a larger scaling 
length.) 



G. Transmitted intensity 

If the characteristics of the initial molecular beam are 
known, the transverse and longitudinal acceptances of 
the full beam transport line may be used to calculate 
the beam intensity at the final focus. As an example, 

which has a skim- 



consider the beamline in section IV F 



mer radius of 0.75 mm and an angular acceptance of 
± 3.3 mrad. From the point of view of an observer at 
the source, the entire beam transport line intercepts the 
same solid angle as would a 2 mm diameter collimator, 
located 0.3 m from the skimmer, with no lenses. The in- 
tensity of an unfocused beam at this location can (often) 
be calculated or measured. From this, and a knowledge 
of the rotational state population fraction, and the ve- 
locity spread, we can estimate the intensity at the final 
focus of the beamline. 

For the methyl fluoride example, we assume a (seeded) 
jet source beam temperature of 3.5 K. (This is based 
upon the equations in Ref. |L8| for a 0.0035 cm diame- 
ter source orifice, a source pressure of 1.9 x 10 5 Pa (1400 
torr) of Ar at a reservoir temperature of 300 K, and a 
methyl fluoride seed of 5%. The source could be either 
pulsed or continuous.) The 3.5 K, results in a J = pop- 
ulation of about 30 %, based upon a Maxwell Boltzmann 
distribution, and a kinetic energy spread of 11.4% which 
implies a longitudinal acceptance of about 60 %. 

Thus an intensity, equivalent to 18 % of a unfocused 
methyl fluoride beam passing through a 2 mm dia. col- 
limator 0.3 m from the skimmer, would reach the fi- 
nal focus some 14.5 m from the source. And, since the 
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FIG. 9: Lens placement of a complete transport line with matching lenses and a 15-m FODO lattice (center line), and the 
calculated rms beam envelope, in x and y, for 560 m/s beam of methyl fluoride in the J = state. We model the beam from a 
skimmed jet source and focus it to a 2 mm diameter spot. A schematic of the beam transport line is shown in Fig. |i[ With a 
central field of 3.2 MV/m, scaling lengths of ro = 6 mm for the FODO lattice lenses, and 12 mm for the matching lenses, the 
beam transport line has an acceptance of 2.5 mm - mrad in both transverse planes. The longitudinal acceptance is about 84% 
for a ± 10 % rms energy spread. 



beam transport line acceptances demonstrate transmis- 
sion without loss, the same intensity should also reach a 
final focus much further from the source. 

Finally, we note that much higher performance beam 
transport lines can be designed. If one uses stronger elec- 
tric fields, shorter FODO cells, fills more of the beamline 
with focusing elements, and increases the size of the final 
focus, most of the solid angle from the skimmer can be 
accepted. Similar measures would allow one to efficiently 
transport and focus faster beams of molecules, such as 
those seeded in a helium jet source. 

V. APPLICATIONS 

Improving alternating gradient transport and focusing 
will make it easier to use molecules in strong-field seek- 
ing states for experiments, for beam transport, or to focus 
molecules for easier detection. In beam resonance exper- 
iments, molecules prepared in a strong-field seeking state 
can be detected, after a transition to a weak-field seeking 
state, as a flop-in resonance. 

The optimized multipole lenses can greatly reduce the 
problem of Majorana transitions (l^, p0| . These are tran- 
sitions that arise because different mj levels belonging 
to the same J are degenerate in zero field. In very weak 
electric fields, a time-varying component caused, for ex- 
ample, by the motion of the molecule through the lens, 
can induce a transition to a different raj state with a 
very different Stark effect. This leads to beam loss, or 
loss of signal and a large background in sensitive reso- 
nance experiments. 

The optimized multipole lenses reduce this problem by 
allowing one to have all normal or all skew F and D lenses 
with their central electric fields in the same direction. 
This is done by changing the sign of 03, the hexapole 
strength (and 05) while leaving the dipole strength (a±) 



unchanged. Field direction changes in going between F 
lenses and D lenses are eliminated. If a weak dipole 
(bias) field is added to the region between the lenses, 
the molecules may never be in a rapidly changing weak 
electric field. 

The problem of Majorana transitions will be greatest 
for molecules in weak-field seeking states focused by pure 
quadrupole and/or sextupole lenses. These lenses have a 
vanishing field at the center and rapid changes in the field 
direction at their entrances and exits 0. The problem 
will be diminished if the molecules are focused in alter- 
nating gradient fields using optimized multipole lenses 
with their central electric fields in the same direction 
and with alternating positive and negative values of 03. 
These lenses have a non-vanishing field everywhere and 
can be optimized for a quadratic Stark effect by choosing 
I a$\ = —0.2. Again a small dipole bias field can be used 
between lenses. 

The problem of Majorana transitions can be eliminated 
by choosing the J = state; which is non-degenerate, 
always strong-field seeking, has the largest Stark effect 
of any rotational level, and is highly populated in a cold 
jet source beam. 

The J = and other strong-field seeking states have 
unique and useful properties that can be exploited in 
experiments. Within a rotational level, J, the \mj\ = J 
states are strong-field seeking. States that are strong- 
field seeking in weak fields remain strong-field seeking in 
stronger fields (but weak-field seeking states will become 
strong- field seeking in the limit of strong fields). Thus, 
for strong field seeking states, there is no restriction on 
the size of the electric field that can be used to focus 
them. This is an advantage for molecules with small 
rotational constants and large dipole moments, where the 
weak field-seeking states in low rotational levels become 
strong-field seeking in modest electric fields. CsF (Fig. [j]) 
and other heavy alkali halides are examples. 
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Long distance beam transport, which can exceed 100 m 
or more, has a number of applications. Since monatomic 
carrier gasses, clusters, and many contaminants will not 
focus through the beam transport line, it can be used 
to clean up a beam. For hazardous and radioactive 
molecules, a long beamline allows one to separate the 
source material and reservoir from the experiment and 
allows one to use radioactive detection in a lower back- 
ground environment. 

The long transit time (54 ms for the 30 m beam line in 
our example) corrected, if necessary for the small longi- 
tudinal velocity changes in the focusing elements, can be 
used for time-of-flight measurements with pulsed sources 
or a beam chopper. The different velocities will focus 
at slightly different longitudinal positions, which may be 
exploited for position-sensitive detection. Alternatively, 
by using a pulsed beam and ramping the electric field in 
the (final focus) lenses, all molecules may be brought to 
a focus, at different times, but at the same position. 

The long flight path may be useful for colincar laser 
excitation of weak transitions. If the colinear laser ex- 
citation is combined with time-of- flight measurement, 



the Doppler spread from the velocity distribution may 
yield information about the absorption profile. Long 
transit time also allows for the decay of some long-lived 
states. Molecular beams may be run in both directions to 
form a very long colliding beam apparatus. In this case, 
the reservoir temperatures of the beam sources may be 
adjusted to equalize the focusing strengths of different 
molecules. 
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